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ABSTRACT 

Hll regions are the birth places of stars, and as such they provide the best measure of current star formation 
rates (SFRs) in galaxies. The close proximity of the Magellanic Clouds allows us to probe the nature of these 
star forming regions at small spatial scales. To study the Hll regions, we compute the bolometric infrared flux, 
or total infrared (TIR), by integrating the flux from 8 to 500 /im. The TIR provides a measure of the obscured 
star formation because the UV photons from hot young stars are absorbed by dust and re-emitted across the 
mid-to-far-infrared (IR) spectrum. We aim to determine the monochromatic IR band that most accurately traces 
the TIR and produces an accurate obscured SFR over large spatial scales. We present the spatial analysis, via 
aperture/annulus photometry, of 16 Large Magellanic Cloud (LMC) and 16 Small Magellanic Cloud (SMC) Hll 
region complexes using the Spitzer Space Telescope 's IRAC (3.6, 4.5, 8 /im) and MIPS (24, 70, 160 /im) bands. 
Ultraviolet rocket data (1500 and 1900 A) and SHASSA Ha data are also included. All data are convolved to 
the MIPS 160 /im resolution (40 arcsec full width at half-maximum), and apertures have a minimum radius of 
35". The IRAC, MIPS, UV, and Ha spatial analysis are compared with the spatial analysis of the TIR. We find 
that nearly all of the LMC and SMC Hll region spectral energy distributions (SEDs) peak around 70 /im at all 
radii, from ~ 10 to ~ 400 pc from the central ionizing sources. As a result, we find the following: the sizes 
of Hll regions as probed by 70 /im is approximately equal to the sizes as probed by TIR (« 70 pc in radius); 
the radial profile of the 70 /im flux, normalized by TIR, is constant at all radii (70 /im ~ 0.45 TIR); the la 
standard deviation of the 70 /im fluxes, normalized by TIR, is a lower fraction of the mean (0.05-0.12 out to 
~ 220 pc) than the normalized 8, 24, and 160 /im normalized fluxes (0. 12-0.52); and these results are the same 
for the LMC and the SMC. From these results, we argue that 70 /im is the most suitable IR band to use as a 
monochromatic obscured star formation indicator because it most accurately reproduces the TIR of Hll regions 
in the LMC and SMC and over large spatial scales. We also explore the general trends of the 8, 24, 70, and 
160 /im bands in the LMC and SMC Hll region SEDs, radial surface brightness profiles, sizes, and normalized 
(by TIR) radial flux profiles. We derive an obscured SFR equation that is modified from the literature to use 
70 /im luminosity, SFR (M yr 1 ) = 9.7(0.7) x 10~ 44 L 70 (ergs s" 1 ), which is applicable from 10 to 300 pc 
distance from the center of an Hll region. We include an analysis of the spatial variations around Hll regions 
between the obscured star formation indicators given by the IR and the unobscured star formation indicators 
given by UV and H<x We compute obscured and unobscured SFRs using equations from the literature and 
examine the spatial variations of the SFRs around Hll regions. 

Subject headings: dust,extinction - galaxies: individual (LMC, SMC) - galaxies: ISM - Hll regions 
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1. INTRODUCTION 

Hll regions are locations of active or recent star forma- 
tion where the extreme ultraviolet (UV) radiation from mas- 
sive OB stars ionizes the surrounding gas (for a good review 
see chapter 5 of Tielens 2005, and references therein). It is 
common to consider Hll regions as being large complexes 
of overlapping individual Hll regions where the sizes and 
components are not determined solely by the classical Strom- 
gren radius of an individual star but by an inclusion of all 
of the physical regimes in the interstellar medium (ISM) af- 
fected by the far-UV photons of the central sources. These 
Hll region components include the central ionizing OB stars, 
and the surrounding photodissociation regions (PDRs) and 
molecular clouds, where the chemistry and t he heating are 
driven by the UV photons from nearby stars dTielensI 12005; 
iRelano & Kennicuttll2009l: IWatson et al.ll2008l) . Thus, Hll re- 
gions in this context are really Hll complexes, or more gener- 
ally, star forming regions. 

The observed properties of these star forming regions, in- 
cluding size and shape, depend on the physical traits such 



as the number of young hot ionizing stars, density of neutral 
gas, abundance of dust, star formation history, and prior su- 
pernoyae (e.g.. see|Hodgdll974tlTielensll2Q Q5: Wal born et all 
12001 ISnideretal.1 1200 9: Harris & Zaritsky 2009). These 
physical traits are often correlated, and the use of data at many 
wavelengths to observe H II regions in different galactic envi- 
ronments is required to piece them together. 

The nature of th e hot central OB stars can be stud- 
ied using UV (e.g., ISmith et alJ Il987t iMartin et al.l 12005 
or nebular em ission lines, suc h as Ha (e.g., lHenizdfl95 
Hod ge & Kennicuttlll983l: iGaustad et al.ll2001l) . The Ha lu 
minosity function of a galaxy is frequently used to determine 
many of the physical properties of Hll regions, including the 
number of ionizing sta rs, evolutionary effect s , and possible 
environmental effects (Ken nicutt et a l. 1989; lOey & C larke 
1998). Hll region size distributions are directly related to the 
luminosity functions and rel ate to the numbers of nebula of 
a give n size for a galaxy (e.g.,[van~den Bergh 1981; Oe y et al.l 
2003). The PDRs can be studied by many methods including 
the [Ol] a nd [Clll infrared (IR) cooling lines and CO radio 
data (e.g., Kaufm an et al.l fl999). The molecular clouds are 
typically studied using radio observations of molecular rota- 
tional lines (e.g., Cohe n et al.lfl98 8: Genzel 1991; Fuku i et al.l 
2008). 

Present star formation rates (SFRs) are calculated using a 
tracer of the UV photons from the y oung mas sive stars and 
spectral synthesis models (see Kennicutt 1998). Hll regions 
will have some fraction of their UV photons obscured by dust 
and some fraction unobscured. For unobscured SFRs, the ion- 
izing photons can be directly observed via UV observations 
or recombination lines such as Ha. For obscured Hll re- 
gions, bolometric IR observations of dust (i.e., the total in- 
frared (TIR)) can be used to recover the extinguished UV 
photons. This is because the dust absorption cross section 
is highly peaked in the UV, and the re-emitted flux is in the 
broad spectral range from the mid-to-far-IR (Kennicutt 1998). 
Because the TIR around Hll regions accounts for all of the ex- 
tincted UV photons, the TIR is expected to be the single best 
indicator of SFR obscured by dust. 

Combining UV, optical, and IR observations of Hll re- 
gions across a whole galaxy allows us to probe the current 
SFR of that galaxy. The Large Magellanic Cloud (LMC) 
and Small Magellanic Cloud (SMC) are ideal natural labo- 
ratories for studying star forming regions and their effects 
on the IS M. The close p roximities of the LMC and SMC, at 
~ 52 kpc (ISzewczvk et a l. 2008) and - 60 kpc dHilditch et al.l 
2005), respectively, allow for detailed star formation stud- 
ies down to parsec or sub-parsec scales depending on wave- 
length. Any broad study of Hll regions across the LMC and 
SMC can take advantage of many multiw avelength observa - 
tions, including the rocket UV data from Smit h et aTl (|1987|) . 
the Southern H a Sky Survey Atlas (SHASSA) Ha data from 
Gaustad et al. (2001), and the Spitzer Space Telescope IR data 
from Meixner et al.l d2006) and K. D. Gordon et al. (2010, in 
preparation). Furthermore, there are many past optical Hll re- 
gion surve ys that c atalog the sources of Ha in the LMC and 
SMC (see iHenizeJ 1956; D avies et al.lll976h iBica & Schmittl 
fl995h . 

Because of the observational advantages, many researchers 
rely on a single-band star formation indicator (i.e., UV, Ha, 
Pa a, 8 /im, 24 /im, etc). There are complications in the UV 
and optical lines in that they can be greatly impacted by ex- 
tinction, thus, requiring extinction corrections. Another com- 
plication is that the observed UV photons can come from 



stars of various ages (< 100 Myr) and will greatly depend 
on galaxy type (i.e., quiescent spiral galaxies, starbursts, etc.; 
(ICalzetti etal.1 12005b ). The 8 and 24 /im IR band emission 
will likely depend on the environment of the host galaxy be- 
cause the abundance of the aromatic s/small grains that give 
rise to their emission depend upon the metallicity and ioniz- 
ing r adiation present (K . D. Gordon et al. 2010, in prepara- 
tion: iDraine & Lil 12007b jDraine et all 120071) . Work done by 
ICalzetti et all (120071) and lDale et al. (2005) indicate that 8 /im 
makes for a poor star formation indicator due to large variabil- 
ity of emission in galaxies with respect to spectral energy dis- 
tribution (SED) shap e and metaUicitv. ICalzetti et all (120071) 
also note, along with ICalzetti et al. (2005), that a star forma- 
tion indi cator using 24 /im by itself can vary from galaxy to 
galaxy. iDale et all (120051) claim that SFRs calculated from 
24 /im emission may be off by a factor of 5 due to variations 
of 24 /im flux with respect to SEDs observed across nearby 
galaxies. 

To compensate for the extinction effects in UV and optical 
nebular emission lines, many researchers are now measuring 
SFRs via a combination of obscured (TIR, 8 /im, 24 /im) and 
unobscured (UV, Ha, Pa a) star formation indicators (e.g., 
Calzetti et al. 2007; Kennicutt et al. 2007; Thilker et al. 2007; 
Relano & Kennicutt 2009; Kennicutt et al. 2009). However, 
the noted differences in 8 and 24 /im emission, relative to 
host galaxy properties, may still introduce uncertainties to 
the calculated SFRs when applying them across large galaxy 
samples. In th eir analysis of 33 galaxies from the Spitzer 
SINGS sample, ICalzetti et all ((2007) claim that a combina- 
tion of Ha and 24 /im gi yes the most robust SFR using 
a procedure similar to the iGordon et al] (l2000|) "flux ratio 
method" 13 . Their calibration has a caveat in that it is useful 
for actively star forming galaxies where the energy ou tput is 
dominated by young stellar populations (Calzetti et al. 2007). 
Kennicutt et al. (2009) analyze SFRs of nearby galaxies de- 
rived by combining Ha with 8 /im, 24 /im, and the TIR. They 
find that linear combinations of Ha and TIR provide for the 
most robust SFRs. 

There is little work done on investigating the effi cacy of 
using 70 or 160 /im as star formation indicators. Dal e et al] 
(2005) claim that the 70 /im emission may make a good 
monochromatic obscured star formation indicator because 
the 70 to 160 /im ratio correlates well with local SFRs. 
A Spitzer analysis of far-IR co mpact sources in the LMC 
(Ivan Loon etal.ll2010al) and SMC (Ivan Loon et al.ll2010bl) . in- 
cluding compact Hll regions, find that the bolometric correc- 
tion to 70 /im is modest, due to a typical dust tempe rature 
of 40 K. In a study of dwarf irregular galaxies, I Walter etail 
(2007) find a good correlation between the brightest 70 /im 
regions and optical tracers of star formation, albeit, with some 
galaxies contributing significant 70 /im diffuse emission at 
large radii. 

Which of the Spitzer IR bands most accurately repro- 
duces the TIR over a large spatial scale? Employing 
aperture/annulus photometry, we analyze 16 LMC and 16 
SMC Hll region complexes using the Spitzer Infrared Array 
Camera (IRAC) and Multiband Imaging Photometer (MIPS) 
bands. We determine that the MIPS 70 /im band provides for 
the most accurate monochromatic obscured star formation in- 
dicator based on an analysis of the Hll region complex SEDs, 
sizes, and radial monochromatic IR fluxes (normalized by the 

13 The "flux ratio method" uses U V and IR fluxes of galaxies to derive 
extinction-corrected UV luminosities (Gordon et al. 2000) 



TIR). We include an analysis of the spatial distribution of the 
unobscured star formation indicators, UV and Ha, relative to 
the IR obscured star formation i ndicators. We m odify an es- 
tablished TIR SFR recipe from Kennicutt (Il998[) to derive a 
new monochromatic obscured SFR equation using the 70 /im 
luminosity. 

In Sectional we list the LMC and SMC Hll region com- 
plexes sampled in this work. In Section [51 we explain the IR, 
UV, and Ha observations and data reduction. The analysis 
of the multi wavelength photometry is discussed in Section H 
The basic results of the photometry, SEDs, and radial pro- 
files, are discussed in Section \5\ as well as a discussion of our 
calculation of the TIR. Discussions of Hll region normalized 
radial SEDs, Hll complex sizes, normalized radial profiles, 
and SFRs are presented in Section [6) In this section we also 
present our derived 70 /im obscured SFR equation. We finish 
with some concluding statements in Section [T] 

2. Hll REGION SAMPLE 

The 16 LMC and 16 SMC Hll regions were selected by 
the following three criteria: the center must be peaked in 24 
/im emission, there must be a nearby peak in Ha, and the 
total sample of Hll regions must sample the full size of the 
LMC and SMC as observed in 24 /im. Along with sampling 
Hll regions across the entire LMC and SMC, the 16 LMC and 
16 SMC Hll regions cover a wide range of sizes and temper- 
atures. The final sample size is a compromise between ob- 
taining good statistics and avoiding overlap between the large 
apertures covering each Hll region. Specifics of the Hll re- 
gions are listed in Table [T] The maximum radii are chosen to 
include the full extent of the Hll region com plex. We describe 
how we quantify this in section I47TI and 16.21 

Shown in Figures [T] and [2] are LMC and SMC images of 
the unobscured and obscured star formation indicators of UV 
(blue-LMC only), Ha (green contours), and 24 /im (red). 
The Hll regions are labeled with circles denoting the max- 
imum extent to which we measure their photometry (see 
Columns 5 and 6 in Table [T]). We do not have UV data for 
the SMC. A mo re de tailed discussion of these figures is pre- 
sented in section \6A\ For full Spitzer IRAC and MIPS images 
of the entire LMC and SMC, see Meixner et al. (2006), K. D. 
Go rdon et al. (2010, in prep aration) , Ivan Loon et al.l (|20 1 0a) , 
and lvan Loon et al.l (l2010bh . 

To view the dust structure of the Hll regions, three-color 
images of the LMC Hll regions are shown in Figure [3] The 
SMC Hll region three-color dust images are shown in Fig- 
ure HI In both figures, the aromatic 8 /im emission is in 
blue, the warm 24 /im dust emission is green, and the colder 
160 /im dust emission is red. The sizes of the Hll regions 
span from 105 pc in radius to 450 pc in radius. Although 
we describe the 8 /im emission as aromatic, we do not at- 
tempt to remove any possible non-aromatic components. The 
IRAC 8 /im band may have some contaminatio n from emis- 
sion of lines such as Hi, fAr iil, and [Arm] (Peeter s~et al.l 
I2002fc iLebouteiller et al.l 12007k . A more detailed discussion 
is presented in section 1531 

3. OBSERVATIONS 

3.1. IR 

The infrared images are created from the Spitzer Space 
Telescope's IRAC and MIPS instruments for the Surveying 
the Agents of a Galaxy's Evolution (SAGE) project (see 
Werner et al. 2004; Meixner et al. 2006; K. D. Gordon et 



TABLE l 
Sample of Hii Regions 



No. Name a 



a (2000) 
(h:m:s) 



5 (2000) 



Maximum Radius 
(") (pc)* 



LMC 



1 


N4 


4:52:08 


-66:55:20 


1120 


280 


2 


Nil 


4:56:48 


-66:24:41 


1505 


380 


3 


N30 


5:13:51 


-67:27:22 


875 


220 


4 


N44 


5:22:12 


-67:58:31 


980 


250 


5 


N48 


5:25:50 


-66:15:03 


1470 


370 


6 


N55 


5:32:33 


-66:27:20 


910 


230 


7 


N59 


5:35:23 


-67:34:46 


875 


220 


8 


N79 


4:51:54 


-69:23:29 


560 


140 


9 


N105 


5:09:52 


-68:52:59 


1785 


450 


10 


N119 


5:18:40 


-69:14:27 


980 


250 


11 


N144 


5:26:47 


-68:48:48 


455 


115 


12 


N157 


5:38:36 


-69:05:33 


1225 


310 


13 


N160 


5:39:44 


-69:38:47 


805 


205 


14 


N180 


5:48:38 


-70:02:04 


1015 


255 


15 


N191 


5:04:39 


-70:54:34 


525 


130 


16 


N206 


5:31:22 


-71:04:10 


1400 


355 


SMC 


1 


DEM74 


0:53:14 


-73:12:18 


350 


105 


2 


N13 


0:45:23 


-73:22:52 


350 


105 


3 


N17 


0:46:42 


-73:31:04 


385 


115 


4 


N19 


0:48:26 


-73:05:59 


385 


115 


5 


N22 


0:48:09 


-73:14:56 


350 


105 


6 


N36 


0:50:31 


-72:52:30 


525 


155 


7 


N50 


0:53:26 


-72:42:56 


560 


165 


8 


N51 


0:52:40 


-73:26:29 


350 


105 


9 


N63 


0:58:17 


-72:38:57 


350 


105 


10 


N66 


0:59:06 


-72:10:44 


700 


205 


11 


N71 


1:00:59 


-71:35:30 


490 


145 


12 


N76 


1:03:43 


-72:03:19 


350 


105 


13 


N78 


1:05:06 


-71:59:36 


350 


105 


14 


N80 


1:08:34 


-71:59:43 


560 


165 


15 


N84 


1:14:05 


-73:17:04 


980 


290 


16 


N90 


1:29:35 


-73:33:44 


385 


115 



Name associated with closest Ha source from Henize i 1956) and 



iDavieset all (19761) . 

Using the measured distances of ~52,000 pc to the LMC 
iSzewczvk et al. 2008) and ~60,500 pc to the SMC (Hilditch et al. 2005). 

al. 2010, in preparation). The observations and data reduc- 
tions are fully descr ibed in the SAGE LMC overview paper 
(Mei xner et al. 2006) and in the SAGE SMC overview paper 
(K. D. Gordon et al. 2010, in preparation). 

IRAC provides MIR imaging data in four pas sbands cen- 
tered around 3.6 /im, 4.5 /im, 5.8 /im, and 8 /im dFazio et al] 
120041) . MIPS provides FIR imaging data in thre e passbands 
cente red around 24 /im, 70 /im, and 160 /im dRieke et al] 
120041) . The LMC and SMC IRAC exposures consist of 
l.°l x l.°l tiles. The LMC MIPS exposures are scan legs 4° 
long, and the SMC MIPS exposures are scan legs 2 ° -5 ° long. 
The exposure depths are the same for both galaxies as a result 
of the same observing strategies, and variations between the 
IRAC tiles and MIPS strips have been removed. The observa- 
tions were taken at two epochs separated by ~ 3 (LMC/IRAC, 
LMC/MIPS, SMC/IRAC) and - 9 (SMC/MIPS) months. 

The full mosaics, for each IRAC and MIPS band, are cre- 
ated via proce ssing using the Wisc onsin pipeline (IRAC; see 
description in Meixner et al.ll2006l) an d the MIPS DAT anal- 
ysis tool (MIPS; Gord on et al.l 12005). The background has 
been subtracted from the full mosaics of each galaxy to re- 
move zodiacal and Milky Way cirrus emission. The mosaics 
cover - 8 ° x 8 ° for the LMC and - 30 deg 2 for the SMC. 

The angular resolution of the IRAC mosaics are l."7, l."7, 
l."9, and 2" for the 3.6, 4.5, 5.8, and 8 /im bands. The an- 
gular resolution of the MIPS mosaics are 6" ', 18", and 40" 
for the 24, 70, and 160 /im bands. So that the data have 




Fig. 1.— LMC image of the UV (blue) from Smith et al. (1987), Ha (green contours) from lGaustad et all <200lh . and 24 fim (red) from Mei xner et al.l 12006). 
UV and 24 jum are convolved to the MIPS 160 /im resolution (40 // ). The SHASSA Ha fluxes are in their original resolution (~ 45"). The selected Hn regions 
are labeled and marked with an aperture that corresponds to the largest spatial extent for which we measure the photometry (see Table [!}• The bar in the lower 
right represents 1 kpc. North is up, and east is to the left. 



equivalent resolution, all of the IRAC and MIPS point-spread 
functions (PSFs) have been transformed to the 40" PSF of 
the M IPS 160 /imb and using the custom convolution kernels 
from lGordonetafl d2008h . 

3.2. ANCILLARY DATA 

The UV data are from imaging tak en with sounding- 
rocket intrumentation (Smi th et all fl987). The bandpasses, 
centered around 1495 A and 1934 A, are 200 A wide full 
width at half -maximum (FWHM ) and 220 A wide FWHM, 
respectively (Smi th et all 11987b . In this paper, we use 
1500 A and 1900 A to refer to these UV bands. The 



1500 A image has a total exposure time of 114s, and the 
1900 A image has a total exposure time of 130s. The 
1500 A PSFs and the 1900 A PSFs have been convolved to 
m atch the 160 /mi PSF using the custom convolution kernels 
of lGordonetal](l2008h . 

The Ha images are from the southern s ky wide-angle imag - 
ing survey known as the SHASSA fsee lGaustad et al.l l2001). 
The SHASSA survey uses a Canon camera with a 52 mm fo- 
cal length lens and an Ha filter centered at 6563 A with a 32 A 
bandwidth. The LMC and SMC images are each comprised of 
five separate 20 m exposures. The LM C and SMC are in fields 
013 and 010 in lGaustadetail (1200 ll) . respectively. SHASSA 




Fig. 2. — SMC image of the Ha (green contours) from Gaustad et al. (2001) and 24 /im (red) from K. D. Gordon et al. (2010, in preparation). The 24 /mi 
fluxes are convolved to the MIPS 160 /mi resolution (40 /; ). The SHASSA Ha fluxes are in their original resolution (~ 45 /; ). The selected Hn regions are 
labeled and marked with an aperture that corresponds to the largest spatial extent for which we measure the photometry (see Table [!}• The bar in the lower right 
represents 1 kpc. North is up, and east is to the left. 



images have a similar resolution (r 
160 /im images (40"). 



'45") to that of the MIPS 



4. ANALYSIS 

4.1. Photometry 

The LMC and SMC photometries are measured for the 
IRAC (3.6 /im, 4.5 /im, 8 /im) bands, MIPS (24 /im, 70 /im, 
160 /im) bands, and the ancillary data for all 32 Hll regions 
spanning physical sizes from — 105 pc in radius to — 450 pc 
in radius (see Tableland Figures [5] &|4]). The photometry of 
the Hll regions is measured in concentric annuli from the cen- 
tral inner 35" aperture outward to the annulus with the largest 
radii. The radii of the annuli are chosen such that they fully 
sample the PSF and are large enough to attain good signal-to- 
noise in the outer parts of the Hll regions where the UV, Ha, 
and 24 /im fluxes quickly drop to background levels. The flux 
for a given annulus is calculated by taking the flux of a larger 
aperture and subtracting off the flux from the adjacent smaller 
aperture. From these annuli, we create radial SEDs and radial 
flux profiles of each Hll region. 

The central apertures for each Hll region are centered 
around the peak 24 /im flux because the peak 24 /im flux 
is observed to closely c oincide with the peak Ha flux 
(Relano & Kennicutt 2009). Also, the 24 /im flux is observed 
to peak around OB stars at spatial scales down to less than 
a parsec (Snid er et~ai] [2009). Each 24 /im peak is visually 



checked to make sure that it coincides with a nearby p eak in 
the Ha S HASSA data and in t he Ha catalog of either lHenizd 
(1956) or Dav ies et all (Il976|) . An Ha peak is considered 
nearby a 24 /im peak if it falls within the central three annuli. 

The 35" apertures correspond to slightly different physi- 
cal sc ales for the LMC an d SMC. The LMC is -52 kpc 
away (ISzewczyk et al.| [2008) which corresponds to a physi- 
cal aperture radius o f w 9 pc. The SMC is —60.5 kpc away 
( Hildi tch et al.l |2005|) which corresponds to a physical aper- 
ture radius of « 10 pc. Thus, we are restricted to measuring 
the properties associated with the dust around Hll regions at 
scales larger than « 9 pc for the LMC and « 10 pc for the 
SMC. 

The largest apertures of the Hll region complexes are de- 
termined, by eye, to be where the coldest gas, emitting at 
160 /im, drops in flux to approximately the level of the back- 
ground diffuse emission or where the Hll region begins to 
overlap another Hll region. The Hll regio n co mplex sizes are 
more quantitatively determined in section [6721 

To remove the sky background and the local LMC or SMC 
background, due to diffuse emission, we perform a back- 
ground flux removal for each Hll region. The background 
flux for every annulus of a given Hll region is computed with 
an annulus taken at radii of 1.1-1.5 times the radius of the 
largest aperture. The flux error for a given annulus is then 
the standard deviation of the background sky flux times the 




Fig. 3.— Three-color dust images of the 16 LMC Hn regions using the IRAC and MIPS data from Meixner et al. (2006). The IRAC 8 /im is blue, the MIPS 
24 /mi is green, and the coldest dust at the MIPS 160 /mi band is red. The 8 and 24 /mi fluxes are convolved to the MIPS 160 /mi resolution (40 /x ). The bars in 
the lower right of each image represent 30 pc. The white circles represent the largest apertures used for the photometry (see Table[T}- North is up and east is to 
the left. The lack of 24 /im flux in the center of 30 Doradus (N157) is due to saturation in the core. 



square root of the number of pixels in the annulus. In the in- 
stances where a neighboring Hll region will fall within the 
background sky annulus, a mask is created that nulls the pixel 
values in the background sky annulus where any large 24 fim 
peaks are observed. The interloping Hll region is nulled in 
every band. The SMC is a smaller galaxy with smaller dis- 
tances between bright Hll regions, relative to the LMC (see 
Figures [Band [2]). This crowding sets up a lower angular size 
limit of the maximum aperture used for many of the SMC Hll 
regions. 

We also perform the total IR photometry of each Hll region 
using an aperture with the largest radius (see Table Q]) to ac- 



quire the cumulative flux. The se fluxes are used to plot the 
Hll region SEDs in section I57T1 For all cumulative fluxes, the 
background flux is taken from the same sky annulus as in the 
photometry of the individual annuli. 

The total IR fluxes of the LMC and SMC galaxies are 
computed using aperture photometry with apertures that en- 
close each galaxy. The LMC circular aperture is centered 

around 5 h :17 m :50 s right ascension and -68°:24 / :42 // decli- 
nation (epoch 2000 coordinates) with a 14,100" radius. The 

SMC aperture is an ellipse centered around 1 .06 m :00 s right 
ascension and -72°:50 / :46 // declination (epoch 2000 coordi- 




Fig. 4. — Three-color dust images of the 16 SMC Hn regions using the IRAC and MIPS data from K. D. Gordon et al. (2010, in preparation). The IRAC 8 /im 
is blue, the MIPS 24 /im is green, and the coldest dust at the MIPS 160 /im is red. The 8 and 24 /j,m fluxes are convolved to the MIPS 160 /im resolution (40" )• 
The bars in the lower right of each image represent 30 pc. The white circles represent the largest apertures used for the photometry (see Table[T}- North is up and 
east is to the left. 

nates) with a major axis radius of 9710" and a minor axis 
radius of 9140". The LMC and SMC sky backgrounds are 
taken from regions outside of their apertures and away from 
any bright IRAC/MIPS infrared sources. 

The number of Hll regions with a given aperture/annulus 
decreases considerably with larger radii for both the LMC 
and SMC. Thus, the conclusions in this work suffer from low 
number statistics at larger radii. For the LMC, the number 
of Hll regions drops to five at an aperture/outer annulus ra- 
dius of 1225". For the SMC, the number drops to five at an 
aperture/outer annulus radius of 525". 



5. RESULTS: 

5.1. SEDs 

Hll region IR SEDs are a combination of the Rayleigh- 
Jeans tail of the stellar photospheric emission and emis- 
sion from dust grain populations. The 3.6 /im fluxes, 
and to a lesser degree the 4.5 /im flu xes, are associate d 
with stellar cont i nua m ore than dust (Hel ouet all 120041) . 
Engelbracht et al. (2008) find that the transition from stel- 
lar dominated emission to dust dominated emission in star- 
burst galaxies occurs around 4.5 /im. The dust emission 
at 5.8 and 8 /im is likely dominated by aromatic emission 



(Puge t & Legeii[T989h . To more conservatively separate the 
dust emission from IR emission associated with stellar pho- 
tospheres, we do not include the 5.8 /im band in our analysis 

14 

The 24, 70, and 160 /im FIR emission are associated with 
silicate and carbonaceous grains of various siz es, as demon- 
strated by interstel l ar dust grain models (e.g ., IDraine & Lee 
[19841; iDesertetafl [19901: iLi & Draind [2001]) . IDraine & Li 
(2007) find that small grains, of sizes ~ 15-40 A, may 
contribute a significant portion of the 24 /im continuum via 
single-photon heating. The 70 /im emission may be due to 
a significant fraction of both small grains and large grains, 
but 160/im fluxes are dominated by large grains (Deser teTaT] 

cm. 

From our annulus photometry, we produce IR SEDs for 
each annulus of each Hn region. Each of these radial Hn 
region IR SEDs is created by using the fluxes (mJy Hz) at 
3.6, 4.5, 8, 24, 70, and 160 /mi. For each Hn region, we an- 
alyze the SEDs and compare how they change radially from 
the core to the outermost measured ann ulus. A full analysis 
of the radial SEDs is discussed in section I67T1 

We also produce IR SEDs using the cumulative flux of each 
Hll region and compare these with the SEDs created using 
the total galactic LMC and SMC fluxes. The total LMC and 
SMC fluxes were comp uted using a single aperture around 
each galaxy (see section PI) . The LMC, SMC, and Hll re- 
gion cumulative fluxes, in Jy, and their la uncertainties are 
tabulated in Table [2J The SEDs are plotted in Figure \5jp for 
the LMC, and Figure [5^ for the SMC. The maximum aper- 
ture radius used for each Hll region correspond to the aper- 
tures marked in Figures [3] and [4] and in the last two columns 
ofTableffl 

The total LMC SED we plot in Figure^ (gray triple-dot- 
dashed line) is very similar, in ter ms of relative strengths be- 
tween bands, to that produced in B ernard et al . (2008). The 
total SMC SED we plot in Figure [5j? (gray triple-dot-dashed 
line) is consistent with the SMC SED in K. D. Gordon et al. 
(2010, in preparation). Nearly all of the LMC/SMC Hll re- 
gions in our sample peak around 70 /im. The few Hll regions 
where the 160 /im flux exceeds that of the 70 /im flux, in 
particular N4, N44, and N180 in the LMC, have maximum 
aperture sizes that are among the largest in the sample. The 
paucity of aromatic emission i n the SMC, relative to the LMC, 
is further discussed in section [63l 

Our results match well with those of Degioia-Eastwood 
(1992), who observed six Hll regions in the LMC using IRAS 
and noted the peak fluxes to be near 60 fim. Other studies have 
found peak wavelengths of Hll regio n SEDs at similar wave- 
lengths, including the results of Indeb etouw et al.1 (120081) . who 
use models of typical Hll regions and find the peak emission 
between 40 and 160 /im, depending on the model used. Simi- 
lar work has b een done on entire gala xies. From modeli ng and 
observations, iDale et al.1 (120011) and iDale et al.1 (120051) show 
that nearby galaxies have IR SEDs that peak between 40 and 
160 /im. 

The large fluxes of the 30 Doradus (N157) SED in Fig- 
ure 02 (red dot-dash line) are due to an abundance of hot 
stars. 30 Doradus is the largest H II region in the Local Group 
(IWalborn et al.ll2002l) at « 200 pc in size (Rubio et al. 1998), 
and c ontains more than 100 OB stars (Walborn & Blades 
1 1997b . Many of the stars are of type 03, which are among 

14 The use of the 5.8 /xm band is further complicated by its relatively low 
sensitivity among the IRAC bands (Fazio et al. 2004; Meixner et al. 2006) 



the h ottest, most luminous stars known dMassey & Hunted 
119981) . 30 Doradus is the brightest Hll region in our sam- 
ple at approximately an order of magnitude fainter than the 
total cumulative LMC flux. There is no comparably bright 
Hll region in our sample, other than perhaps N160, which is 
spatially located near the 30 Do radus comp lex. An analysis 
of 30 Doradus in the FIR by lAguirre et afl (120031) finds that 
this Hll region contains ~ 20% the total LMC FIR emission. 
The brightness we compute is a lower limit because both the 
24 /im and 160 /im fluxes are saturated in the core. The sat- 
uration in the central annuli at 24 /im artificially lowers the 
surface brightness we measure out to ~ 35 pc. The saturation 
in the central aperture of the 160 /im emission lowers the sur- 
face brightness we measure out to ~ 9 pc. We do not include 
30 Doradus in our quantitative analysis of star formation in- 
dicators, but we do incl ude a brief qualitative discussion of 
this Hll region in section [63] in order to highlight some of the 
effects a hotter Hll region may have on our conclusions. 

The SMC Hll region IR SED fluxes are about an order 
of magnitude, or more, fainter than the LMC Hll region IR 
SED fluxes. There are several possible explanations for this. 
For example, metal poor galaxies will have relatively less 
dust extinction giving rise to a greater interstellar radiation 
field and less IR emission, or an abundance of very small 
grains/aro matics can preferentially increase the Wien' s side of 
the SED (iGalliano et al] 12005 iGalametz et al J 120091) . How- 
ever, there are systematically fewer ionizing photons (Ha) 
and fewer reprocessed IR photons (TIR) in the S MC Hll re- 
gions relative to the LMC Hll regions (see Section [6.4. II and 
Figure [TIT) TIR and Ha are good star formation indicators 
(Kennicuttl ll998l) . and decrements in both of them are evi- 
dence that intrinsically lower SFRs are occurring in the SMC 
Hll regions. 

This can also be understood by considering the luminos- 
ity functions for Hll regions. The Ha luminosity function 
follows a power law where the number of Hll regions with 
a given range of luminosities is proportional to L~ a (e.g., see 
IKennicutt & Hodgdl98fl[l986l:IOev & CJa^e1IT998l) . In gen- 
eral for galaxies, the shape of the Ha luminosity function de- 
pends on the number of ionizing photons and the effects from 
evolution (lOev & Clarkel 119981) . IKennicutt & Hodgd fl986) 
find an a - 1.65 ±0.15 for the LMC and an a - 1.75 ±0.15 
for the SMC using - 20 LMC and ~ 100 SMC Hll re- 
gions JromJhe_£ata]ogof JD^IeietaD (Il976|) . Recent work 
from lGielesetal.1 (120061) finds that the LMC and SMC have 
power laws of a ~ 2. There is no universal upper cutoff for 
Hll region Ha luminosity functions, so galaxies with higher 
SFRs will statistically have more luminous Hll regions. Be- 
tween the LMC and the SMC, IKennicutt etaD (11989b find 
more of the most lumino us Hll regions in the LMC, and 
IKennicutt & Hodgel (119861) explicitly show the LMC Ha lu- 
minosity function cuts off at larger luminosities, which is not 
surprising considering the presence of 30 Doradus. There 
are few calculations of Hll region IR luminosity functions, 
so it is uncertain how the Ha luminosity function compares 
with the IR luminosi ty function. However, in one such study, 
iLivanou et al.1 (l2007|) compare the frequency distribution of 
IRAS 100 /im emission for LMC and SMC Hll regions, and 
find that the LMC Hll regions have a higher IR luminosity 
cutoff (see their Figure 4). 

5.2. TIR 

The bolometric IR flux, TIR, is a measure of the total 
dust emission in the IR. TIR makes for a sensitive tracer of 
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Fig. 5 . — IR SEDs of the cumulative fluxes (mJy Hz) of the Hn regions (calculated using the largest aperture) in the -a) LMC and -b) SMC. The total galactic 
LMC and SMC fluxes are also plotted as triple-dot-dashed lines in gray. The computation of the Rayleigh- Jeans tail (160-500 fim) is explained in Section 5.2. 



obscured star formation around Hll regions becau se the ab- 
sorpti on cross section of dust is highest in the UV ( Ken nicuttl 
119981) . We calculate TIR for each Hll region annulus, the cu- 
mulative Hll region apertures, and the total LMC and SMC 
galaxy apertures. Our calculation of TIR includes the 8, 24, 
70, and 160 /im fluxes, as well as the Rayleigh- Jeans tail, 
from 160 to 500 /im, of a computed modified blackbody 
(MBB) function. We do not use the 3.6 or 4.5 /im bands so as 
to avoid contributions from stellar continua. 

The TIR calculation is done by integrating the area under 
the linear interpolation of the flux (from 8 to 160 /im), in 
flux (F v ) versus frequency (y) space, then adding the far- 
IR Rayleigh- Jeans tail of the MBB. The total flux under 
the Rayleigh- Jeans tail is computed by integrating under the 
MBB, in flux (F v ) versus frequency (v) space, from 160 to 
500 /im and using frequency steps equivalent to 0.25 A. 

The MBB is defined as 



MBB = X ^ B U (T) 



(1) 



where x is a constant, /3 is the emissivity index, and B is the 
Planck function (mJy) with dust color temperature T (K). The 
shape of the MBB component is calculated using the ratio of 
the 70 and 160 /im fluxes to numerically derive a tempera- 
ture of the dust 15 , an d a = 2. We choose the classical emis- 
sivity index of /? = 2 dGezari et al. II 19731: iDraine & Leelll984l) 
for simplicity. The emissivity index depends upon the chem- 
ical nature and internal physic s of the grains and is generally 
thought to be from 1 < /3 < 2 (Dupa ceFaT]l2003h with some 
indications that it can b e greater than two (Mennella et al.1 
119981: iDesert et all 120081) . Determining beta is problematic 
because T and beta are degenerate. Also, flux errors can make 
determining beta more problematic (Shetty et al. 2009). 

15 The larger grains associated with the 70 and 160 /im fluxes are least 
affected by single-photon heating and, thus, produce a more accurate Planck 
function. 



The uncertainty in TIR is calculated via standard error prop- 
agation and using the uncertainties in the 8, 24, 70, and 
160 /im bands as well as the uncertainty in the MBB. The 
uncertainty in the MBB component is calculated via a Monte 
Carlo method. Random 70 and 160 /im fluxes are generated 
assuming Gaussian distributions about the measured fluxes 
with an FWHM w 235cr F . The new MBB dust color tem- 
perature is then calculated using the ratio of the new 70 and 
160 /im Planck functions. The flux from 160 to 500 /im is 
integrated for this new MBB component as described before. 
This process is repeated 1000 times, and the standard devia- 
tion of the integrated fluxes is the la error in our MBB mea- 
surement. The MBB uncertainty is assumed to be indepen- 
dent of the other flux errors despite being dependent upon the 
70 and 160 /im fluxes. Thus, our final quoted TIR uncertain- 
ties are slightly larger than the true TIR uncertainties. There 
are additional uncertainties not taken into account in our cal- 
culations such as differences in the emissivity index. We also 
do not know the contribution of single-photon heating to the 
70 or 160 /im fluxes that might lead us to calculate a temper- 
ature that is too high. 

The MBB component (160-500 /im) is not the dominant 
contributor to the IR bolometric flux because al l of the SED s 
are peaked blueward of 160 /im (see sections l5Jl and [67Tb . 
The hottest Hll region is 30 Dor (N157) at 29.2 K, which puts 
the peak of the MBB at ~ 100 /im. The coolest Hll region is 
N180 at 22.4 K, which corresponds to a peak in the MBB at 
~ 130 /im. The LMC Hll regions are generally cooler than 
the SMC Hll regions, using the cumulative fluxes from Ta- 
ble [2] However, the LMC cumulative fluxes are measured 
using larger apertures which introduce a greater percentage of 
flux from cold dust at greater radii. 

The individual annuli of each Hll region exhibit a wider 
range of temperatures than the cumulative fluxes in Table [2] 
The LMC Hll region annulus with the warmest dust is around 
30 Dor, with a color temperature of 41.5 K and a peak in the 
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TABLE 2 














IR Fluxes of Hii Regions 








Name 


3.6/im 


4.5/xm 


8/xm 


24/xm 


70/im 


160/im 


■ l c 




(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(K) 


LMC 


LMC^ 


2150 ±40 


1280 ±20 


6860 ±140 


8080 ±320 


127000 ±6400 


265000 ±32000 


23.9 


N4 


13.4 ±0.24 


9.88±0.19 


57.3 ±1.2 


49.5 ±2.0 


742 ±37 


2040 ±240 


22.6 


Nil 


17.8 ±0.33 


13.6±0.3 


198 ±4 


238 ±10 


3730 ±190 


7700 ±920 


23.9 


N30 


4.44 ±0.10 


3.43 ±0.08 


28.2 ±0.6 


27.2±1.1 


483 ±24 


1100 ±130 


23.4 


N44 


10.2±0.19 


7.71 ±0.15 


133 ±3 


204 ±8 


2930 ±150 


4420 ±530 


25.6 


N48 


9.32 ±0.23 


6.65 ±0.16 


101 ±2 


78.9 ±3.2 


1670 ±80 


4560 ±550 


22.6 


N55 


4.48 ±0.10 


3.27 ±0.07 


16.2 ±0.3 


22.2 ±0.9 


377 ±19 


629 ±75 


25.1 


N59 


6.56±0.15 


5.09±0.11 


63.0±1.3 


101 ±4 


1490 ±70 


2570 ±310 


24.9 


N79 


3.65 ±0.07 


2.99 ±0.06 


40.4 ±0.8 


81.6±3.3 


1060 ±50 


1290 ±160 


26.9 


N105 


42.8 ±0.8 


27.0 ±0.5 


259 ±5 


213±9 


4620 ±230 


8360 ±1000 


24.6 


N119 


9.59 ±0.23 


7.22 ±0.17 


120 ±3 


120 ±5 


2510±130 


3530 ±420 


26.0 


N144 


2.19 ±0.04 


1.73 ±0.03 


23.7 ±0.5 


29.4 ±1.2 


345 ±17 


800 ±96 


23.4 


N151 b 


46.9 ±0.9 


45.9 ±0.9 


550 ±12 


2340 ±90 


22100 ±1100 


19500 ±2300 


29.2 


N160 


21.2±0.4 


17.3 ±0.3 


277 ±6 


702 ±28 


8820 ±440 


8910 ±1070 


28.2 


N180 


7.13 ±0.21 


5.16±0.16 


60.0 ±1.3 


44.3 ±1.8 


947 ±47 


2690 ±320 


22.4 


N191 


1.23 ±0.03 


0.892 ±0.021 


7.22±0.15 


14.1 ±0.6 


152 ±8 


235 ±28 


25.5 


N206 


3.78±0.14 


4.28 ±0.11 


92.2 ±1.9 


77.1 ±3.1 


1550 ±80 


3300 ±400 


23.8 


SMC 


SMC 


330 ±10 


240 ±10 


240 ±10 


350 ±10 


10800 ±500 


20000 ±2400 


24.5 


DEM74 


0.521 ±0.027 


0.355 ±0.017 


1.31 ±0.03 


2.24 ±0.09 


77.4 ±3.9 


152±18 


24.2 


N13 


0.474 ±0.046 


0.420 ±0.034 


3.45 ±0.09 


7.16±0.30 


154 ±8 


270 ±33 


24.8 


N17 


0.876 ±0.030 


0.640 ±0.020 


1.63 ±0.04 


2.64 ±0.11 


61.3±3.2 


123 ±15 


24.1 


N19 


1.47 ±0.05 


1.09 ±0.03 


4.54±0.10 


14.1 ±0.6 


356 ±18 


502 ±60 


26.0 


N22 


1.13 ±0.03 


0.848 ±0.024 


4.33±0.10 


18.2 ±0.7 


255 ±13 


343 ±41 


26.3 


N36 


0.937 ±0.043 


0.780 ±0.031 


2.53 ±0.06 


11.8±0.5 


321 ±16 


419 ±50 


26.5 


N50 


2.17 ±0.05 


1.67 ±0.04 


2.13 ±0.05 


5.66 ±0.23 


163 ±8 


200 ±24 


26.9 


N51 


0.0960 ±0.0240 


0.0570 ±0.0150 


0.779 ±0.019 


2.08 ±0.08 


57.3 ±2.9 


116±14 


24.0 


N63 


0.545 ±0.020 


0.384 ±0.013 


0.558 ±0.013 


1.56 ±0.06 


34.0 ±1.7 


50.9 ±6.1 


25.7 


N66 


4.13 ±0.09 


3.29 ±0.08 


5.34±0.11 


38.9 ±1.6 


481 ±24 


531 ±64 


27.6 


N71 


1.39 ±0.03 


1.02 ±0.02 


1.67 ±0.04 


2.29 ±0.09 


56.1 ±2.8 


128 ±15 


23.4 


N76 


1.49 ±0.05 


1.10±0.03 


2.17 ±0.05 


6.91 ±0.28 


116±6 


204 ±25 


24.8 


N78 


1.20 ±0.03 


0.945 ±0.024 


2.05 ±0.04 


14.2 ±0.6 


105 ±5 


136 ±16 


26.5 


N80 


0.836 ±0.027 


0.579 ±0.019 


1.21 ±0.03 


2.46 ±0.10 


95.6±4.8 


191 ±23 


24.1 


N84 


2.29 ±0.06 


1.77 ±0.04 


6.48 ±0.14 


20.1 ±0.8 


413 ±21 


835 ±100 


24.0 


N90 


0.226 ±0.005 


0.177 ±0.004 


0.311 ±0.007 


0.847 ±0.034 


20.0 ±1.0 


38.9±4.7 


24.2 



Note. — Fluxes and derived dust color temperatures for the largest aperture of each Hii region (see Column 5 in Table0. ±lcr uncertainties 
include the statistical and calibration errors summed in quadrature. The IRAC 3.6, 4.5, and 8 /im calibrations have 1.8%, 1.9% and 2.1% uncertain- 
ties, respectively (Reach et al. 2005). The MIPS 24 /im calibration has a 4% uncertainty lEngelbracht et al. 2007), the MIPS 70 /im calibration has 
a 5% uncertainty I Gordon et al. 2007), and the MIPS 160 /im calibration has a 12% uncertainty I Stansberrv et al. 2007). 
a Dust color temperature calculated using the 70 and 160 /im fluxes (see section |5^2l . 

30 Doradus 24 /im and 160 /im fluxes are saturated in the core. 



MBB at ~ 70 /im. The LMC Hii region annulus with the 
coldest dust is around N180, with a color temperature of 21.3 
K and a peak in the MBB at - 135 /im. The SMC Hii region 
annulus with the warmest dust is around N36, with a color 
temperature of 39.2 K and a peak in the MBB at ~ 75 /mi. 
The SMC Hii region annulus with the coldest dust is around 
N90, with a color temperature of 19.1 K and a peak in the 
MBB at ~ 150 /im. The annuli with the warmest dust are 
always the central aperture, except for N36, N50, and N63 in 
the SMC. These Hii regions have the warmest dust in annuli 
just outside of the center aperture. The coldest dust for each 
Hii region is always found at the outer annuli. The SED s of 
the Hii region annuli are further discussed in Section |6J"1 

For the LMC Hii regions, the MBB integrated component 
ranges from ~ 2% of the 8-160 /im integrated fluxes for the 
inner annuli of 30 Doradus to ~ 25% for the outer annuli of 
the cooler Hii regions (N48, N180, and N206). For the SMC 
Hii regions, the percentage of the MBB integrated fluxes to 
that of the 8-160 /im integrated fluxes range from ~ 4% (in- 
ner aperture of N66) to as high as ~ 30% for the coldest Hii 
region outer annulus (N90). 

If the peak of a blackbody or MBB shifts slightly redward 
or blueward, the intensity near the peak will not significantly 



change. However, even slight shifts in the peak of the black- 
body will greatly affect the intensities on the far Wien or 
Rayleigh- Jeans sides of the blackbody. As noted above, all 
of the Hii region SEDs in our sampl e pea k between ~ 70 and 
150 /im, and as discussed in Section l63l the 70 /im emission 
is « 40% -50% of the TIR, on average, for both the LMC and 
SMC at all radii measured. Therefore, the MBB Rayleigh- 
Jeans component we measure is greatly dependent upon the 
70 /im flux. The spreads in the peak of the SEDs give rise to a 
spread in the relative importance of the Rayleigh- Jeans tail to 
the TIR. The fluxes in the Rayleigh- Jeans tails of the MBBs, 
relative to the TIR, range from ~ 1.7%, for a dust color tem- 
perature of 41.5 K, in the central annuli of 30 Doradus, to 
~ 25%, for a dust color temperature of 19.1 K in the outer 
annuli of N90. 

Our defined TIR (8-500 /im) is similar to other calcula- 
tions of IR bolometric flux in the literature, except we do not 
g o quite as far into the far-IR as some (e.g., 8-1000 /im used 
in lKennicuttl [l998). Most of the IR flux around obs cured star 
formi ng regions will emit in the 10-120 /im range (Ken nicuttl 
119981) . We believe our TIR results are robust; we do not as- 
sume a single SED parameterization but instead compute the 
TIR for each individual Hii region annulus. Our TIR values 
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are ne arly identical to what are predicted by iDale & H elou 
(120021) for mo s t of th e annuli around our Hll regions. The 
IDale & Heloul (120021) formalism (see Equation (4) in their 
work) is derived using IRAS and Infrared Space Observa- 
tory (ISO) observations of normal galaxies and using the 60 
to 100 /im flux ratios to parameterize the range of possible 
galaxy SEDs. Our results differ by a factor of a few in the 
outer annuli of our Hll regions. This is not unexpected be- 
cause the cold dust at 160 /im contributes more at these radii. 
IDale & Heloul (2002) make a point of noting that their TIR 
results will be less accurate for colder galaxies because the 
parameterization is calculated using the IRAS 100 /im band 
which does not extend to long enough wavelengths to include 
the flux from the coldest dust. 

5.3. Morphology and Radial Surface Brightness Profiles 

The morphology of the dust around the Hll regions is 
displayed in Figures [3] and [4] There is much in the litera- 
ture pertaining to H l l region gas and dust mor p hologies (see 
Relano & Kennicutt 2009; Snider etal. 2009: IWatson et al.1 
20081: [Smith & Brooks 2007); iHelou et alJf2004 . and our re- 
sults are similar. The warm dust at 24 /im is heavily peaked 
around the centers of the Hll regions and near peaks in Ha 
emission (see Figures [T] and [2]). The 24 /im emission drops 
off quickly, whereas, the aromatic 8 /im emission is more 
extended. The cold dust measured by 160 //m is very ex- 
tended and dominates the outer environs of the Hll region 
complexes. On smaller spatial scales than is resolved in this 
work, others have noted that 8 /im is largely absent near the 
hot OB stars in the center of Hll regions, possibly due to the 
destructio n of aromatics or low dust density near the inner 
OB s tars (Snider et al. 2009; Wa tson et al.ll2008l: IHelou et al.l 
120041) . iLebouteillereta l. (2007), by observing NGC 3603, 
provide spectroscopic evidence that aromatics are destroyed 
near the central cluster, whereas, very small grains survive 
nearer the central ionizing source. 

Within the limiting resolution of this work, the dust emis- 
sion at every wavelength peaks near the centers of the Hll 
regions, as we show by plotting the LMC and SMC Hll re- 
gion surface brightness radial profiles in Fig ure [6 ] The TIR 
is the bolometric IR flux calculated in section [5^21 and the la 
uncertainties are on the order of the sizes of the points. The 
scatter can be understood as being due to intrinsic differences 
in the SEDs of the Hll regions. 30 Doradus (N157) has the 
highest surface brightness profile at all wavelengths, as shown 
in Figure^ (boxes). The surface brightnesses radial profiles 
plotted here wi ll be used to examine the broad Hll complex 
sizes in section lo\2l 

The single-band IR surface brightness profiles in Figure [6] 
are similar to the bolometric (TIR) surface brightness profiles 
in that they all peak in the center of the Hll regions, even for 
the coldest dust at 160 /im. The central surface brightness is 
lower for the SMC Hll regions than the LMC Hll regions by 
about an order of magnitude, for all bands. This is possibly a 
resolution effect. The resolution of our work, using our small- 
est aperture, is ~ 20 pc in diameter. Any Hll region smaller 
than this size, i.e., not resolved, will have an artificially low 
surface brightness. This might conceivably affect both the 
LMC and SMC, but the SMC Hll luminosity function has a 
slightl y steeper slope and a lower high mass cutoff than the 
LMC (IKennicutt & Hodgdl 19861) . Thus, the SMC may have 
a larger number of Hll regions that are small enough to be 
unresolved. To better study this, a larger sample of LMC and 
SMC Hll region surface brightnesses, of similar luminosities, 



will need to be analyzed. Also, we need a better understand- 
ing of how IR Hll region luminosity functions correlate with 
the Ha luminosity functions. 

The 8 /im fluxes peak in the center despite studies show- 
ing a lack of aromatics near the hot central OB stars and 
an abund ance of aromatics on the surface of outlying PDRs 
(Rela no & Kennicuttl 120091: ISnider et al.l 120091: IWatson et al.l 
2008). Our inner aperture resolution of ~ 20 pc diameter is 
not as high as many of those previous studies. We likely are 
not able to resolve the destruction or removal of the aromatics 
near the central hot stars, and geometry may play a role via 
foreground emission. There is also the possible contamina- 
tion in the 8 /i m band of non- aromatics, such as Hi, J Aril"] , 
and [Arm] fsee lPeeters et al.ll2002l: iLebouteiller et al.ll2007l) . 
It is unlikely that these other lines dominate the emission at 
8 /im, but they may contribute a significant fraction. This 
may be particularly true in the SMC where there is a paucity 
of aromatics. 

6. DISCUSSION 

6. 1 . Normalized Radial SEDs 

From the results of section 15.11 we find that most of the 
Hll regions in the LMC and SMC have IR SEDs that peak 
around 70 /im. Thus, we expect 70 /im to be a dominant con- 
tributor to the bolometric IR flux. However, there is intrinsic 
scat ter i n the radial profiles of all of the IR bands (see sec- 
tion [53]). To better understand the spatial variations of the IR 
band emission in LMC and SMC Hll regions, we analyze the 
radial SEDs at radii from near their core (~ 10 pc) to their 
outskirts. Because we are interested in how the IR dust bands 
relate to the bolometric flux, we use the radial SED fluxes nor- 
malized by the TIR, where the normalized flux at frequency v 
is defined as 

F N ,„ = z//v/TIR. (2) 

In Figure [71 we show the average normalized (by TIR) flux, 
in units of F N ,„, for each annulus in the LMC (panel a) and 
SMC (panel b). Each line represents an average of all of the 
Hll regions with a given annulus. For both the LMC and 
SMC, the inner annuli show relatively strong average F N ^4 
values that drop as the annuli proceed radially outward. The 
average F N ,7o values in both the LMC and SMC remain rel- 
atively constant, with slight variations. The average F N ,i6o 
values for both the LMC and SMC behave opposite to that of 
the average F NJ 24 values in that they increases radially outward 
from the center of the Hll region. 

The normalized 8 /im aromatic emission behaves very dif- 
ferently in the LMC than in the SMC. In the LMC the av- 
erage F N ,8 values increase slightly with radius. In the SMC 
the average F N ,8 values do not follow an obvious trend with 
radius. The 8 /im emission in the SMC Hll regions on av- 
erage accounts for about 7% of their TIR, while in the LMC 
Hll regions, aromatics account for 15% of the TIR. The lower 
normalized 8 /im values in the SMC are at least partially due 
to the lower intrinsic fluxes at 8 /im for SMC Hll regions, 
and the SMC in general (see Figure [5]). The lack of a radial 
trend for the 8 /im emission in the SMC may be a combi- 
nation of the lack of aromatics, in which we would expect 
a trend similar to the LMC 8 /im where aromatics emission 
increases away from the central ionizing source, and poten- 
tial contributions to the 8 /im from other bands. For instance, 
the [Aril] nebular emission line is not expected to contribute 
much to the 8 /im emission (Leboute iller et al .1120071) . but the 
lack of strong aromatics in the SMC means that this weaker 
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Fig. 6. — Radial surface brightness profiles of the -a) 16 LMC Hn regions and the -b) 16 SMC Hn regions. The surface brightnesses are in units of 
10 30 [ergs s _1 pc -2 ], and the aperture radii are in units of [arcseconds] with spatial units of [parsecs] on the top axis. The panels, from top to bottom, are TIR, 
8 /im, 24 /im, 70 fim, and 160 jum. The ±lcr statistical errors are on the order of the point sizes. The b oxes in panel (a) are the points associated with the 
measured surface brightness profile of 30 Doradus (N157). The original LMC MIPS 70 /im data are from Meixner et al. (2006), and the original SMC MIPS 
70 jum data are from K. D. Gordon et al. (2010, in preparation). 
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FIG. 7. — Average normalized radial IR SEDs for the -a) LMC and -b) SMC. Each line represents the average of the normalized (by TIR) fluxes of all of the 
Hn regions in an annulus with inner and outer radii (labeled in the upper left). The parenthetical numbers to the left of the labels correspond to the number of 
Hn regions included for each annulus. The colors of the lines go from the inner annuli (blue), to middle annuli (green and red), and finally to the outer annuli 
(black and gray). 
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line, along with Hi and [Arm] emission (Pee ters et al.ll2002t 
Lebouteiller et al. 2007), may be more significant factors in 
the observed 8 /im fluxes. The contribution of non-aromatic 
emission to the 8 /im fluxes within SMC Hll regions needs 
to be more fully explored. The lower SMC 8 /im emission is 
further discussed in section [631 

The F N ,7o values are nearly constant among all annuli for 
both the LMC and SMC (70 /im « 45% TIR). Thus, the SEDs 
in each annuli must all peak near 70 \im and the peak of the 
SED must not considerably vary as we spatially observe the 
Hll regions from their centers outward. We further explore 
the implications of this result in the following sections. 

6.2. Physical Sizes of Hll Region Complexes 

The appropriate size aperture to use for extragalactic 
studies of Hll regions is unknown and is typically set by the 
limiting resolution of the instrument because many of the Hll 
regions are not resolved in the IR at extragalactic distances. 
An incorrect aperture size will add uncertainties to any SFR 
calculation because the accumulated flux may not all be from 
the Hll region of interest. We explore the Hll complex size 
scales, as probed by the 8, 24, 70, and 160 /im dust emission, 
and compare them to the size scales found using the TIR flux. 

We quantitatively probe the sizes of the LMC and SMC 
Hll regions by tracing the TIR, 8, 24, 70, and 160 /im ra- 
dial profiles (see Figure [6]) and calculate the radius at which 
each band's surface brightness drops by 95% (~ 2a) from that 
measured in the inner 35" (~ 10 pc) radius aperture 16 . All 
four dust bands peak in the center aperture for every Hll re- 
gion (see Figure [6]) , and the central aperture is where we are 
most certain that young stars dominate the UV output. 

The results of our size determinations can be found in Ta- 
ble [3] The size for each LMC and SMC Hll region complex 
is listed along with the radius in arcseconds where the surface 
brightness drops to less than 5% of the surface brightness of 
the inner aperture. The mean, standard deviation, and frac- 
tional error (standard deviation divided by the mean) for the 
LMC and SMC Hll complexes are also tabulated. The Hll 
region sizes have large standard deviations about the mean 
values. The errors, as a fraction of the mean, are 31% -39% 
for all bands. 

For the TIR in the LMC Hll regions, the mini- 
mum/maximum radii measured are 140/385 arcsec (35/97 pc), 
respectively. For the TIR in the SMC Hll regions, the mini- 
mum/maximum radii measured are 140/525 arcsec (41/154 
pc), respectively. The computed sizes listed in Table [3] are not 
drastically different for the LMC and the SMC, particularly 
given the large standard deviations. 

The sizes calculated from TIR, though not identical, are 
close to the sizes calculated using the 70 /im emission. This 
can most clearly be seen in Figured] where we plot the Hll re- 
gion 8, 24, 70, and 160 /im sizes, normalized by the TIR sizes. 
The mean F N ,7o size falls nearly on the vertical dotted line as- 
sociated with TIR for both the LMC and SMC. The standard 
deviation in the sizes measured using F N ^o is smaller than for 
the sizes measured using the other normalized bands. This is 
a another conseq uence of 70 /im being near the peak of dust 
SED (see Section [67Tb . The mean size of the Hll regions us- 

16 We choose to use a common cutoff for all of the bands. A 3a cutoff 
was attempted, however, it was too strict for many of the Hll regions. The 
surface brightnesses rarely dropped to below 3cr of the central aperture for 
the coldest dust at 160 /xm. This is due to the relatively shallow slopes of the 
radial 160 /xm surface brightnesses 



TABLE 3 
Sizes of Hii Region Complexes 





TIR 


8/im 


24/im 


70/im 


160/im 


Name 




Radius (arcsec) 




LMC 


N4 


140 


210 


105 


140 


210 


Nil 


245 


490 


210 


245 


420 


N30 


210 


280 


140 


210 


350 


N44 


315 


385 


210 


315 


350 


N48 


315 


385 


140 


315 


420 


N55 


245 


245 


210 


245 


245 


N59 


175 


245 


140 


175 


245 


N79 


140 


175 


105 


140 


210 


N105 


175 


210 


140 


175 


245 


N119 


385 


490 


280 


420 


490 


N144 


210 


280 


175 


210 


315 


N157 


NA a 


280 


NA a 


315 


NA a 


N160 


175 


245 


140 


210 


490 


N180 


315 


490 


245 


315 


490 


N191 


140 


140 


105 


140 


140 


N206 


315 


350 


210 


315 


385 


Mean 


233 


306 


170 


243 


334 


Std dev 


79 


113 


54 


82 


114 


Fr.error c 


0.34 


0.37 


0.32 


0.34 


0.34 


SMC 


DEM74 


NA^ 


NA b 


245 


NA^ 


NA b 


N13 


245 


210 


140 


245 


350 


N17 


280 


280 


210 


280 


315 


N19 


280 


280 


245 


280 


315 


N22 


245 


245 


140 


245 


315 


N36 


455 


385 


385 


455 


490 


N50 


525 


490 


315 


525 


560 


N51 


280 


245 


140 


280 


350 


N63 


210 


245 


210 


210 


210 


N66 


210 


210 


175 


245 


280 


N71 


140 


280 


140 


140 


210 


N76 


315 


NA^ 


245 


315 


NA^ 


N78 


140 


140 


105 


140 


210 


N80 


350 


280 


210 


315 


385 


N84 


315 


280 


210 


315 


350 


N90 


175 


210 


140 


175 


210 


Mean 


278 


270 


203 


278 


325 


Std dev 


107 


84 


74 


105 


105 


Fr.error c 


0.39 


0.31 


0.36 


0.38 


0.32 



Note. — Hii co mple x sizes are computed using the method 
described in section fO] 

a Not calculated due to flux saturation in the core of the Hii 
region. 

Not calculated due to a shallow radial profile; the maximum 
Hii region aperture is too small. 
c Fr.error = Std dev/Mean. 

ing either TIR or 70 /im is 60 ± 20 pc the LMC and 80 ± 30 
pc for the SMC. 

The calculated 8 /im sizes in the LMC are as large as or 
larger than sizes measured using 70 /im or TIR, with the ex- 
ception of 30 Doradus (N157). There is no such trend in the 
SMC where the sizes traced by the 8 /im emission may be 
larger or smaller than the sizes calculated using the 70 /im 
dust emission. This is shown in Figure [5] where the sizes of 
the Hii regions, as measured by F N ,g, are systematically larger 
than TIR in the LMC but roughly scattered around TIR in the 
SMC. This trend is consistent with the 8 /im flux following a 
shallower decline with radius in the LMC relative to the SMC. 
The mean size of the Hll regions using 8 /im is 75 ± 35 pc for 
the LMC and 80 ± 25 pc for the SMC. 

All of the Hii regions are smallest when calculated using 
the 24 /im emission. This is consistent with others who have 
found the warm small dust grains probed by 24 /im being 
more heavily peaked near the center of star forming regions 
(e.g., Relano & Kennicutt 2009). None of the Hii regions 
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Fig. 8. — LMC and SMC Hn region sizes of the 8, 24, 70, and 160 /im fluxes normalized by the sizes computed using TIR. The LMC is plotted in the left panel, 
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deviations. All of the Hn regions in Tableware included, but many of the individual Hn regions overlap due to discrete 35 " increasing radii for successive 
annuli. The horizontal dotted line marks the positions where the sizes of the monochromatic bands are equal to the computed sizes measured using TIR. 



have sizes measured via F N ,24 in Figure [5] that are larger than 
the sizes calculated via TIR. The mean size of the H II regions 
using 24 /am is 40 ± 15 pc for the LMC and 60 ± 20 pc for the 
SMC. 

The mean size of the Hll regions using 160 /im is 85 ± 30 
pc for the LMC and 95 ± 30 pc for the SMC. We have shown 
in Figure that the F N ,i6o averages in the Hll re gion SEDs 
are largest at large radii (also discussed in section [63]) . Con- 
sequently, the Hll region sizes, as probed by 160 /im emis- 
sion, are generally larger than the sizes calculated using the 
other IR dust bands. When the sizes measured via the TIR 
and 160 /im emission are compared (see Figure [5]), the sizes 
calculated using the cold 160 /im emitting dust are system- 
atically larger or the same size as the TIR in both the LMC 
and SMC. The sizes measured by the cold dust emission at 
160 /im are the largest despite the fact that the 160 /im sur- 
face brightnesses peak in the center of the Hll regions (see 
Figure©. The increasing importance of the 160 /im emission 
to the TIR in the outer radii of the Hll regions causes a rela- 
tively shallow decline in flux with radius. The computed large 
Hll region cold dust sizes are a consequence of the pervasive- 
ness of the cold dust emission to large radii. 

From this analysis, we conclude that the sizes of Hll re- 
gions, as probed by dust, depend greatly on the wavelength 
observed but also depend greatly on the individual Hll region 



due to the large deviations in each band. The warmer dust 
probed by 24 /im will give smaller Hll region sizes than the 
cold dust probed by 160 /im. The sizes probed by 70 /im 
emission is nearly identical to the sizes probed by the TIR, 
with little scatter. This is a consequence of the IR SEDs peak- 
ing near 70 /im at all s patial distances for nearly every Hll 
region (see section [67Tb . Thus, when taking into account the 
limits of resolution and aperture size for extragalactic stud- 
ies, 70 /im is the most ideal of the monochromatic IR star 
formation indicators because it will trace the TIR spatially 
and extend to physically larger distances than, say, the 24 /im 
emission. Furthermore, we show in Figure^ that for most of 
the Hll regions in our sample, an aperture radius of ~ 100 pc 
would be adequate for acquiring most of the 70 /im flux. 

The size distribution of the nebular portion of Hll regions, 
as traced by Ha, has been extensively studied and argued to 
follow either a power law (Kennicutt & Hodge 1980) or an 
exponential law (Ivan den Bergh 1981). Given that the lumi- 
nosity function is calculated to follow a power law, Oey et al. 
(2003) argue that the size distribution must also follow a 
power law because the luminosity measurements come from 
ionizing photons withi n a specific nebu l ar vo lume. For the 
Magellanic Clouds, iKennicutt & Hodgel ([1986) use an expo- 
nential function to estimate the scale length for Hll region 
sizes to be ~ 80 pc. From our data, we cannot determine a 
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statistically significant size distribution for the LMC or SMC 
because our sample is too small. The initial selection of our 
32 Hll regions was by eye which will bias any size distribu- 
tions created from this data. A more objective IR selection 
of a larger number of LMC and SMC Hll regions is required 
for an analysis of the LMC and SMC Hll region size distri- 
butions, particularly at the small end. Ignoring the warmer 
dust probed via 24 / mi, the sizes we do measure are roughly 
the same as that of the Kennicutt & Hodge (1986) derived Ha 
size distribution scale length. 

It must be stated that the sizes we measure include the gen- 
eral star forming regions. The dust emission is high out to 
relatively large distances. Our size estimates are not directly 
comparable to the Ha-derived nebular sizes. Although there 
may be similarities in that the IR-measured sizes of the LMC 
Hll regions should be somewhat larger than those in the SMC 
because the Ha and IRAS 100 /xm luminosity functions in the 
LMC have a hi gher mass cut-off than in the SM C (see discus- 
sion i n Section [5TT1 Kennic utt & Hod ge 1986: iLivanou et all 
120071) . The large standard deviations of the IR-derived sizes 
is real because the scatter in the surface brightn ess p rofiles 
from which they are derived is real (see section 15.31) . This 
scatter is probably introduced by sampling Hll regions of dif- 
ferent high mass stellar populations and evolutionary histo- 
ries, whi ch are two com ponents that shape the Ha luminosity 
function (lOev & Clarkell 19981) . 

6.3. TIR Normalized Radial Profiles 

Plotted in Figure[9]are the 8, 24, 70, and 160 /xm radial 
profiles, normalized by TIR, (F N ,„), of the LMC and SMC 
Hll regions. For the LMC Hll regions, the F N ,8 emission is 
truncated near the center of the Hll regions (down to ~ 10% 
TIR) and is relatively constant (~ 20% -25% TIR) for annuli 
with radii greater than ~ 200" (~ 50 pc). The F NJ 24 emission 
is sharply peaked near the centers of the Hll regions (~ 20% 
TIR) and falls to a modest level (~ 5% - 10% TIR) for annuli 
with radii greater than ~ 200" (~ 50 pc). The F N ,7o emission 
is relatively constant for all radii probed here (~ 40% - 50% 
TIR). The F N ,i6o emission increases dramatically with aper- 
ture/annuli radii, from (~ 20% TIR) near the core to (~ 50% 
TIR) near the outer annuli. The scatter about the mean for 
F N ,i6o is larger than for the other emission bands. The 8 /xm 
outliers in Figure ^jfi are from the N4 Hll region, where two 
bright 8 /xm sources are included in the outer annuli at ~ 260 
pc from the center (see N4 in Figure [3]). 

For the SMC Hll regions, the F N ,8 does not show an in- 
crease from the central core as it does in the LMC. The F N ,8 is 
consistently around 5% TIR at all radii, which is significantly 
smaller than in the LMC. The F N ,24 emission is peaked in the 
center as it is for the LMC, although to a smaller magnitude 
(~ 15% TIR). The F N ,7o emission in the SMC is consistent 
with that for the LMC, at all annuli radii. The F N ,i60 emission 
is also similar to that for the LMC. 

Due to the change in the shape of the F NJ 24 profile at low 
radii, 24 /xm may not be as robust a tracer of TIR (or obscured 
SFR) as 70 /xm, which has a relatively constant profile with 
radius. An observer would need to develop a conversion from 
24 /im to TIR that is dependent on the radius of the aperture 
used for the measurement. 

The average SMC F N ,8 profile in Figure [9] is depressed rel- 
ative to the LMC at all scales. The same is true, although to a 
lesser extent, for the F N ,24 profile. The 8 and 24 /xm fluxes are 
at least partially associated with aromatics/small grains under- 
going stochastic heating (Puget&Leger 1989; Draine & Li 



l2007h . The low values of F N ,8 and F N ^4 in the SMC rela- 
tive to the LMC are an indication that the SMC has fewer of 
these aromatics/small grains. This is consistent with the re- 
sults in K. D. Gordon et al. (2010, in preparation) who find 
that the SMC has an aromatic fraction of 1.1%, which puts 
the aromatic abund ance in the SMC on the low end of the 
iDraine et afl (|2007) Spitzer SINGS sample. This is also con- 
sistent wi th the SMC's l ow metallicity and high radiation field 
hardness (iGordon et alJl2008h . The strong dependence of the 
host galaxy's aromatic/small grain population presents a prob- 
lem for using either of these bands as a monochromatic star 
formation indicator. Previous works have noted some of these 
potential probl ems when using 8 or 24 /xm as star formati on 
indicators (e.g. JDale et aT1l200llCaIzetti et al.H200ll2007h . 

In Figure [ToJ we show the fractional error (la/mean) of 
the F n , v measurements for each annulus. Plotted are the frac- 
tional errors with respect to annuli radius for the LMC (panel 
a) and SMC (panel b). For all radii, except for perhaps at the 
largest radii measured here where low number statistics dom- 
inate, the fractional error of the F N ,7o indicator is lower than 
the other IR dust emission bands and is ~ 5% - 12% the mean 
value at radii less than 240 pc in the LMC and SMC. The frac- 
tional error of F N ,24 is > 0. 15 until the largest radii where low 
number statistics become an issue. 

Using the IR fluxes, normalized by the TIR, we find that the 
70 /xm emission is a robust tracer of the TIR, and thus the star 
formation obscured by dust. This is because the normalized 
70 /xm flux is relatively constant (70 /xm ~ 0.40-0.50 TIR) 
for all radii measured, does not systematically vary between 
the LMC and SMC, and has the smallest fractional errors at 
nearly every radii. These results are consequences of the fact 
that the SEDs of the H II region annuli in our sample all peak 
near 70 /xm (see section 1670) . 

These results will be less robust for galaxies that have Hll 
regions that are hotter than the typical LMC or SMC Hll 
region. However, we do not expect this to be the case be- 
cause the hottest Hll region in our sample, 30 Doradus, would 
fall on the h igh end of any Hll re gion Ha-derived luminos- 
ity function (Kenn icutt et aljfl989b . If hotter Hll regions are 
more common in a particular galaxy, they will have dust pro- 
files more like 30 Doradus (see boxes in Figure [9)2) where 
the 24 /xm flux is a larger percentage of the TIR out to larger 
radii. Consequently, the normalized 70 /xm flux is reduced 
near the core. This produces an increasing slope in F N ,7o until 
some radius where the 24 /xm emission is no longer a domi- 
nant contributor to the TIR. It is probable that we would see 
this effect on smaller scales than is probed in this work for all 
Hll regions. 

The photons of hotter Hll regions will also affect the ob- 
served 8 and 160 /xm normalized fluxes. The hotter Hll re- 
gions may push the photoionization front to larger radii caus- 
ing a decrement of aromatics near the core. Thus, the peak 
normalized 8 /xm flux would occur at larger radii. The 160 /xm 
flux would not dominate the TIR until even larger radii than 
is measured in this work. If hotter Hll regions are more com- 
mon in other galaxies, then there will be a larger scatter about 
the average F N ,7o ratio with larger variations near the core. 

6.4. Star Formation Rates 
6.4. 1 . Radial SFR Densities 

For the star formation obscured by dust, we use th e 
TIR SFR computation from Equation (4) of lKennicuttl (119981) . 
The [Kennicutt] (11998b obscured SFR equation is valid for 
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starbursts, but dust emission in Hll regions and starburst 
galaxies of varyin g metallicities and ionization measures be - 
have similarly dGordon et al.M2QQ8b lEngelbracht et al.H2008l) . 
Dust obscured star formation in Hll regions is similar to 
star formation in starbursts because starbursts are usually 
very dusty and, furthermore, we re likely the birthplaces 
of most of the pre sent-day star s (lElbaz & C esarskv 2003; 
iBernard-Salas et al.ll2009h . The iKennicuttl (119981) obscured 
SFR equation takes the form of 



SFR = 4.5 x 10" 44 L 



(3) 



where SFR is the star formation rate in units of solar masses 
per year. The value 4.5 x 10~ 44 is a constant derived from 
spectral synthesis models, assumptions on the shape of the 
initial mass function (IMF), and assum ptions on the timescale 
of the star formation event (|Kennicutt 1998). L TIR is the TIR 
luminosity in units of ergs per secon d. Similarly, we use 
Equations 1 and 2 in IKennicuttl (|1998|) to calculate the unob- 
scured SFRs derived from UV and Ha, which rely on similar 
assumptions. In all three cases a Salpeter IMF is assumed 
(Sal D^etel[T955h 

It is noted by IRelano & Kennicuttl (12009b and iHelou et al.l 
(2004) that Ha emission peaks near the maximum 24 /im 
emission. However, the peak UV emission is offset from 
Ha and 24 /im, presumably because the UV emission is ab- 
sorbed by the dust where th e H a and 24 /im emit the strongest 
(Relano & Kennicutt 2009). We also find that Ha peaks near 
the peak 24 /im emission (see Figures[T]and[2]). The UV fluxes 
in the LMC data indicate a peak UV emission near the peak 
Ha and 24 /im fluxes, but as a result of our spatial resolu- 
tion we cannot determine if the UV is offset from the peak 
Ha or 24 /im emission as is claimed by IRelano & Kennicuttl 
(2009). There are regions in the LMC with strong UV fluxes 
but weak or absent Ha and 24 /im fluxes. Many of these UV 
fluxes are found in the centers of noted Hi supergiant shells 
where hot OB stars and supernovae remnants have cleared out 
much of the gas and dust a nd triggered star for mation on their 
peripheries fsee lKim et aBl999UCohen et al.ll2003b . 

In Figure \TT\ we show the radial dependence of SFR den- 
sities (M yr _1 pc -2 ) derived using the obscured SFR indi- 
cator (TIR) and the unobscured SFR indicator (Ha) for the 
LMC (panel a) and the SMC (panel b). The magnitudes of 
the SFR densities depend on the aforementioned assumptions 
about the IMF, star for mation time scales, and spectral synthe- 
sis models used (Kennicutt 1998). However, the magnitudes 
are not crucial for this analysis. It is only the dependence 
with radius that we wish to convey. The general trends of the 
obscured and unobscured SFRs with radius will be constant 
because all conversions from a star formation indicator to a 
SFR depend upon a constant value. 

The obscured and unobscured SFR densities in the LMC 
are larger than in the SMC, which is expect ed g iven the larger 
LMC Hll region luminosities (see Section [57Tb and the high- 
end cutoff at larger luminosities for the LMC Ha and IR lu- 
minosity functions (Kennicutt & Hodge 1986; Livanou et al.l 
l2007h . 

For the obscured SFR densities, both the LMC and SMC 
show a peak in the inner annuli and a decrease outward. This 
is expected given that the TIR emission peaks in the centers 
of the Hll regions as plotted in Figure |6] The unobscured 
SFR densities, calculated via Ha, also peak in the center for 
LMC Hll regions. Because the 24 /im emission peaks where 
the TIR emission peaks in Figure [fO these results are further 



evidence, along with the results in Figure \T\ that the 24 /im 
emission traces the Ha emission. Along the same lines, the 
8, 70, and 160 /im peak fluxes will also trace the peak Ha 
fluxes. These results should break down on smaller scales 
than measured here because the dust grains may be destroyed 
when near the central ionizing stars. For example, it is noted 
that the aromatics probed with 8 /im are not observed to be 
any n earer the central OB stars than in the outlying PDRs 
(Relano & Kennicutt 20091 ISnider et al.l 120091: IWatson et al.l 
2008; Helou et al. 2004). 

The Ha-derived unobscured SFR densities are slightly off- 
set from the TIR-derived obscured SFR densities in the SMC 
(see middle panel of Figure [TT]). There is a slight decline 
of the Ha determined unobscured SFR densities in the inner 
aperture. This may well be a consequence of our limited Ha 
resolution. We aim to explore this in the future with higher 
resolution data. 

As shown in the bottom panel of Figure [Tib , the ratios 
of obscured SFR densities to the unobscured SFR densities 
slightly increase with radius in the LMC. If this trend is real 
then it indicates that Ha is more strongly peaked in the center 
of an Hll region than TIR. One possible explanation is that 
the TIR is a measurement that includes dust at many temper- 
atures. As shown earlier, the colder dust dominates at larger 
radii. This trend is not seen in the SMC, although, we do not 
trace many of the SMC Hll regions out to large radii. 

We do not have UV data of the SMC, but we do have UV 
data of the LMC. For the LMC, we can do the same analy- 
sis of unobscured SFRs via UV as we did with Ha (see Fig- 
ure [12]) • In general, the UV-derived SFR densities are similar 
to the TIR-derived SFR densities in that they are both cen- 
trally peaked. The spatial congruence is probably a resolu- 
tio n effect because it is contr ary to the higher resolution work 
of (Relano & Kennicutt 2009) who found UV emission offset 
from IR emission. There is no discernable trend in the ratio 
of obscured to unobscured SFRs in the bottom panel. 

6.4.2. Obscured SFR Indicators 

It i s well argued tha t the TIR is th e best obscured SFR in- 
dicator (iKennicuttll 19981) . and, indeed. [Kennicutt et al.l (2009) 
find that TIR is a better star formation indicator than 8 or 24 
/im. However, if an observer does not have multiple IR bands 
then it may be difficult to make an accurate calculation of TIR. 
Thus, a monochromatic star formation indicator that can reli- 
ably recover the TIR over a large spatial scale is valuable. 

We substitute the TIR luminosity in Equation [3] with the 8, 
24, 70, and 160 /im luminosities by using the average lumi- 
nosity of a given band, normalized by TIR. The resulting SFR 
equation contains observable quantities. These equations take 
the same form as Eq. [3l but there are band specific constants 
added (Cg, C24, C70, and C160) that are from the factor L v Lj\ R 
used to transform L TIR to L v . From these constants we com- 
pute the SFRs (M yr _1 ) using the measured 8, 24, 70, or 
160 /im luminosities (ergs s _1 ), 



SFR = 4.5xlO- 44 C, y L„ 



(4) 

44 is 



where SFR is in units o f solar mass per year, 4.5 x 10 
the constant from Kennicutt (1998), C v is the IR band specific 
constant, and L v is the observed luminosity in ergs per second. 
In Table ffl we show the aperture size-dependent constants 
used for the aperture radial sizes from 10 to 300 pc for the 
LMC (top rows) and from 10 to 200 pc for the SMC (bottom 
rows). Listed in the columns are the size of the aperture radius 
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TABLE 4 
SFR Equation Average Constants 



Radius 
(pc a ) 



C 8 



C24 



C70 



Cm 



LMC 



10 


8.84±1.80 


5.07 ±0.80 


2.19±0.15 


6.11 ±2.40 


15 


20 


8.30±1.79 


5.89±0.91 


2.10±0.11 


5.42 ±1.70 


15 


50 


6.65 ±1.22 


7.92 ±1.43 


2.08±0.12 


4.05 ±0.93 


15 


100 


5.79 ±0.97 


9.85 ±1.98 


2.11±0.15 


3.39 ±0.74 


15 


200 


5.09 ±0.83 


12.46 ±2.46 


2.18±0.17 


2.83 ±0.65 


13 


300 


4.57 ±0.27 


14.36±2.16 


2.23±0.11 


2.49 ±0.32 


4 


SMC 


10 


16.64 ±3.14 


8.72 ±3.33 


2.22±0.18 


4.46 ±1.62 


16 


20 


16.45 ±3.26 


10.16±3.45 


2.18±0.15 


4.09 ±1.32 


16 


50 


16.70 ±4.40 


13.55 ±4.75 


2.17±0.15 


3.53 ±0.87 


16 


100 


16.49 ±5.45 


16.33 ±5.63 


2.19±0.15 


3.20 ±0.63 


16 


200 


18.82±3.18 


ll.93zb4.61 


2.18±0.19 


3.49 ±1.05 


2 



of 8 and 24jim between the LMC and SMC, as discussed 
in section 16.31 can be seen in the differences of Cg and C24 
in Table 01 The 70 /im band specific constants, C70, do not 
vary between the LMC and SMC, or with radius, within the 
±lcr standard deviations. Any SFRs derived from the 8, 24, 
or 160 am luminositi es should be used with caution. 

The ICalzetti et al.l (|2007|) SFR formalism is commonly 
used in the literature (e.g., IWu et al.l 120081: [Galamet z et al I 
2009|). Using 2 4 /im , we compare the SFRs calculated from 
Kennicutt (1 1998b modified method with the extinction- 



Note. — Constants to the SFR equations are calculated for each aperture radius 
(parsecs). The top values are for the LMC and the bottom values are for the SMC. The 
constants are defined in Section [6.4.21 30 Doradus (N157) has been removed from this 
analysis due to saturation in the core in 24 and 160 fim. The uncertainties are given as 

±1(7. 

a Using the measured dist ances of ^52,000 pc to the LMC I Szewczvk et al. 200|) and 
~60,500 pc to the SMC iHilditch et alJ2005l) . 

in parsecs, the SFR equation constants, and the number of 
Hll regions included in the calculation. 30 Doradus was not 
included in these calculations due to the saturation of 24 and 
160 /im. 

The SFR equations for each of the observable bands are de- 
pendent on the physical scale of the instrumental PSF, partic- 
ularly for the 8, 24, and 160 /im luminosities. The deviation 



our 

corrected method of ICalzetti et al.l (I2007L see their Equation 
(9)) for apert ures of increasing size. The SFR equation in 
ICalzetti et al.l (|2007|) is computed by finding a best-fit line 
through a plot of Pa a luminosity versus 24 /im luminosity, but 
using the highest metallicity galaxies in their sample. In the 
inner radi i (^10 pc) of our sa mple, the results using the for- 
malism of ICalzetti et al.l (|2007|) produce SFRs that are greater 
by a factor of 2 in the LMC and 1.7 in the SMC. The SFRs 
between the two methods slowly equilibrate with radius until 
they are nearly identical at a radius of ~ 280 pc in the LMC 
and ~ 140 pc in the SMC. The results di verge again in the out- 
skirts of the Hll regions. At 200 pc, the ICalzetti et al. I (120071) 
method produces SFRs that are ~ 2/3 the value of our SFRs 
in the LMC and - 3/4 the value of our SFRs in the SMC. It 
is reasonable to expect our SFRs to be lower near the center 
of Hll regions because we do not take into account the un- 
extincted UV photons that is in the ICalzetti et al.l (120071) for- 
malism. The L MC and SMC differe nces are likely because 
the formalism of Calze ttTet al.l d200 7) assumes high metallic- 
ity systems, whereas, the LMC and SMC have different, and 
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low, m etallicities. Thus, the SFR formalism of ICalzetti et alJ 
(2007) is not as accurate for lower metallicity systems, but 
our derived SFR formalism is less accurate if there is a large 
percentage of unextincted UV photons. Thus, we stress that 
o ur star formation formalism provides obscured SFRs. 

IWu et al.l (120081) also use 24 jam as a star formation indi- 
cator. They do a study of 28 compact dwarf gal axies and 
calcul ate SFRs using the 24 /xm SFR recipe from IWu et al.l 
(2005), which was derived using the observed correlations 
between Hq, 1.4 G Hz, and 2 4 /im l uminosities. Unlike the 
ICalzetti et all (120071) sample, IWu et al.1 (2005) include both 
metal-poor an d metal-rich galaxies. The SFRs used via the 
IWu et al.l (120051) method better ma tch the SFRs estimated us- 
ing 1.4 GHz, and lWuetaD (120081) note that this method de- 
rives SFRs that are up to several times la rger than the SFRs 
derived us ing the formalism from lCalzetti et al.l d2007|) . Com- 
paring our lKennicuttl (119981) modified 24 /im me thod of calcu- 
lating radia l cumulative SFR s with the results of ICalzetti et al.l 
(120071) and lWu et al.l (120051) . we find that the SFRs measured 
using our metho d are ~ 35% - 57% the values found using the 
IWu et al.l J2005I) m ethod, for all radii. This differs from the 
results of ICalzetti et al.l (120071) . where our SFRs became rela- 
tively greater at larger radii. Thus, SFRs calculated at larger 
radii using our formalism lie b etween those giving higher 
SFRs using the lWu et all (120051) method and the lower SFRs 
calculated using the ICalzetti et al.l (|2007|) method. These re- 
sults are consistent because the LMC and SMC have lower 



metalli cities than the normal star forming galaxies used to de- 
rive the Calzet tTet al.l ([2007) formalism. 

By observing 22 starburst galaxies. [Brandl et ail (2006) ex- 
plore the use of aromatics as star formation indicators. They 
find a correlation between the aromatic 6.2 /im emission and 
the TIR, where the TIR calculated is uncertain to within a fac- 
tor of 2. Using our measured aromatic 8 /im emission and the 
uncertainties in Cg in Table 01 we can derive TIRs that are un- 
certain by a factor of - 1.2- 1.5 in the LMC and ~ 1.2-2.0 
in the SMC. The differences in our uncertainties are due to the 
changing dispersions of Cg between the LMC and SMC and 
between radii. The SMC correlations are more uncertain in 
our sample than the LMC correlatio ns. Our results are com- 
parable or better than the results of Bran dl et al.l (|2006|) . Our 
measurements do not suffer as much from the introduction 
of diffuse dust emission because we center our obse rvations 
around star forming areas. The Bra ndl et ail (2006) dataset 
uses the Spitzer IRS slits, and the smallest linear scales they 
probe are ab out 200 pc, whereas, we probe down to w 20 pc. 
iBrandl et"aTI ([2006) note that any local variations may be "av- 
eraged out" in their measureme nts. 

The results of Bran dl et all ([2006) illustrate some of the 
difficulties in estimating TIR associated with star formation 
when averaging over an entire galaxy. Our data are acquired 
by observing around the central ionizing sources of the Hll 
regions. Thus, a caveat is that our results are more applicable 
to star forming regions and not entire galaxies, where the IR 
emission of cold diffuse dust is more prominent. As an ex- 
ample, the average C70 value for the LMC and SMC is 2.17, 
which corresponds to L TIR = 0.46 x L70. The derived values of 
C70 using the total LMC and SMC galactic apertures are 2.56 
and 2.22, respectively. The SMC value is on the high end 
measured in Table [H and the total LMC C70 constant is larger, 
beyond the dzlcr standard deviations, of those measured for 
the average Hll regions at any radii in Table H The uncer- 
tainty on the final calculated SFR after integrating over an 
entire galaxy will depend upon how much of the IR photons 
originate from star forming regions and how much are due to 
photons from the general stellar population. 

6.4.3. The 70 \im Star Formation Proxy 

The results in Table [H indicate that 70 /im is the most 
suitable obscured SFR indicator. The C70 constant remains 
the same value for all radii and between the LMC and SMC, 
within the la standard deviations. For these reasons, we 
calculate one mean value of C70 using all of the C70 values 
at every radius from 10 to 300 pc, including both the LMC 
and SMC. From this, we have computed a monochromatic 70 
/im SFR e quation using the formalism from Equation (4) of 
Kennicutt (1998) and the mean value of C70. The monochro- 
matic obscured SFR equation is 



SFR = 9.7(0.7) xlO~ 44 L 



70, 



(5) 



where SFR is in units of solar masses per year and L70 is in 
units of ergs per second. The la standard deviation in Equa- 
tion (0 is from the uncertainties in the spread of C70. The 
mean of the C70 radial averages is 2.17 and the 1 a stan- 
dard deviation (0.15) comes from the mean of the C70 radial 
standard deviations. The la standard deviation of C70 is only 
~ 7% of the mean value of C70. As is the case f or the origi- 
nal obscured SFR equation from Kennicutt] (Il998|) (see Eq.[5]), 
this equation works as a good approximation for the total SFR 
if the area measured is almost entirely obscured by dust. 
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We calculate the total obscured SFR of the LMC, using 
Equation \5\ and the cumulative 70 /im flux of the LMC (see 
section O). The LMC obscured SFR of 0.17 (Afo yr" 1 ) 
agrees well with the results of Harris & Zaritsky (2009) who 
estimate the total mean LMC SFR to be 0.2 (M© yr" 1 ), with 
variations over the last five billion years of up to a factor of 
two. 

The usefulness of Equation © extends to Hll regions mea- 
sured with aperture radii anywhere from 10 to 300 pc. Care 
should be taken if measuring sizes outside of these ranges. On 
very small physical scales our relation, as all others, should 
break down because the effe cts of individual protosta rs and 
stars will become important (Indebetou w et al.ll2008(> . Fur- 
thermore, the effects of hotter Hll regions, like 30 Doradus, 
may result in a radially dependent 70/TIR ratio that would 
introduce errors into our formalism at lower radii (see the 
discussion in section [63]) . There is also the possibility that 
the luminosities of smaller Hll regions are different than ex- 
pected from a presumed IMF due to small number statistics 
(lOev & Clarke 1998). For the same reason, Indeb etouw et al.l 
(2008) argue that measurements of many smaller clusters will 
produce SFRs that are too low since they are not fully sam- 
pling the stellar IMF. 

Equation Q also may not hold at very large radial distances 
because of dust emission excited from older (colder) stellar 
populations contaminating the IR SEDs. Further work must 
be done to model the contribution of older stellar populations 
to the 70 /xm emission, particularly in the outer regions. A 
careful background subtraction will help mitigate any con- 
tamination, and a modest annulus of radius ~ 100 pc is all 
that is required to obtain most of the bolometric dust emis- 
sion in the LMC and SMC Hll regions (see section lo\2l) . Fur- 
thermore, our equation is valid over the wide range of Hll 
region luminosities and sizes probed in this work. The low 
end of the Magellanic Cloud luminosity functions includes 
clusters with luminosities on the order of the Orion Nebula 
(|Kennicutt & Hodgelll986|) . approximately equivalent to N90 
in the SMC, while the high end covers the massive 30 Do- 
radus Hll region. However, this research would benefit from 
a sample of Hll regions in a larger galaxy sample so as to 
more broadly explore environmental effects. 

7. CONCLUSIONS 

We performed UV, H a, and IR aperture/annulus photom- 
etry of 16 LMC and 16 SMC Hll regions at the common 
resolution of the MIPS 160 /xm 40" FWHM PSF. With the 
aperture photometry, we computed IR SEDs using the total 
cumulative Spitzer IRAC (3.6, 4.5, 8 /xm) and MIPS (24, 70, 
160 /xm) fluxes of the LMC, SMC, and each Hll region. The 
cumulative SED of the LMC is brighter than the SMC, and 
their shapes are consistent with IR SEDs of the Magellanic 
Clouds in the literature. The SMC Hll regions have cumula- 
tive fluxes that are weaker, by about an order of magnitude, 
than the LMC Hll regions. We argue that this is expected 
from what is known about the LMC and SMC H II region lu- 
minosity functions. 

We compute the IR bolometric flux (TIR) for each Hll re- 
gion aperture/annulus by integrating under the linear extrap- 
olation of the 8, 24, 70, and 160 /xm fluxes. We include the 
input from the Rayleigh- Jeans tail of a computed dust MBB 
by integrating in steps of 0.25 A, from 160 to 500 /xm. 

Using annulus photometry of the Hll regions, at distances 
of ~ 10 up to ~ 400 pc from their cores, we plot the radial 
8 /xm, 24 /xm, 70 /xm, 160 /xm and TIR surface brightness 



profiles of the 32 Hll regions. The surface brightness profiles 
peak near the center of the Hll regions for all IR bands and 
TIR. We argue that resolution effects may play a role in reduc- 
ing the SMC Hll region surface brightness profiles relative to 
the LMC Hll region surface brightness profiles. The surface 
brightness profiles are used to compute the sizes of each Hll 
region, using the TIR, 8 /xm, 24 /xm, 70 /xm, and 160 /xm 
emission. The typical sizes measured have large scatter and 
are « 75 pc for 8 /xm, « 50 pc for 24 /xm, « 70 pc for both 
70 /xm and TIR, and « 90 pc for 160 /xm. The large scatter for 
each size measurement likely stems from sampling different 
parts of the LMC and SMC luminosity functions. 

We compute normalized (by TIR) IR SEDs for each Hll re- 
gion annulus. The radial SEDs nearly all peak around 70 /xm 
at all radii for every Hll region, out to ~ 400 pc. The result of 
this is that it gives four favorable reasons to choose 70 /xm as 
a monochromatic obscured star formation indicator, and they 
are (1) 70 /xm emission most closely traces the size of Hll re- 
gions as found using the TIR (« 70 pc in radius); (2) 70 /xm 
flux, normalized by TIR, remains nearly constant with radius 
(L 70 ~ 45% L TIR ) from - 10 to 400 pc; (3) 70 /xm flux, nor- 
malized by TIR, has the smallest fractional error (0.05-0.12 
out to 220 pc); and (4) 70 /xm flux, normalized by TIR, does 
not systematically differ between the LMC and SMC. 

Radial SFR density plots, using the obscured indicator 
(TIR) and unobscured indicators (UV, H ap, are computed us- 
ing the SFR recipes of lKennicuttl (I1998|) . The computed SFRs 
for all indicators peak near the center of each Hll region. The 
SFRs for the LMC are larger than for the SMC due to the 
larger luminosities of the LMC Hll regions. Within the lim- 
its of our resolution, the spatial similarities between UV, Ha, 
and 24 /xm are similar to results in the literature. 

A modified v ersion of the obscured SFR equation from 
iKennicuttl (Il998|) is created that depends on a conversion fac- 
tor from TIR to a monochromatic IR band. The conversion 
factors for the 8, 24, and 160 /xm bands all have large un- 
certainties and are dependent on radius and/or host galaxy. 
The conversion factor for 70 /xm is constant, within the low 
uncertainties, at distances from 10 to 300 pc, and between 
the LMC and SMC. We produce a final modified obscured 
SFR equation using a single averaged 70 /xm conversion con- 
stant. This K ennicuttl (1998) modified equation is applica- 
ble for Magellanic-like Hll regions and for aperture sizes of 
10- 300 pc radius. 

The LMC and SMC have different environmental properties 
than those galaxies of earlier Hubble types. The dust emission 
observed around Magellanic Cloud Hll regions may not be 
indicative of dust emission around Hll regions in other galax- 
ies. The analysis in this work would be greatly benefited by 
including Hll regions in a larger sample of galaxies. 

By analyzing the 30 Doradus IR properties, we found pos- 
sible caveats in using 70 /xm as an obscured star formation 
indicator. If Hll regions in a particular galaxy have, on aver- 
age, similar fluxes as 30 Doradus, then the 70 /xm emission, 
relative to the TIR, will exhibit more scatter. However, we do 
not expect this to be the case because 30 Doradus is on the 
extreme h igh end o f Hll region Ha-derived luminosity func- 
tions (iKennicutt et al.lfl9~8 9). 

Another general caveat to mention is that this work applies 
directly to Hll regions and their immediate surroundings. It is 
general practice for those studying SFRs over entire galaxies 
to use the SFR recipes calculated from Hll regions, despite 
the added uncertainties of going from a local star forming re- 
gion to an entire galaxy. The results from this work, applied 
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to entire galaxies, will be subject to similar uncertainties. 

The recent launch of the HERSCHEL telescope will greatly 
increase our understanding of star formation in galaxies. At 
3.5 m, not only is HERSCHEL the largest space telescope ever 
launched, but the PACS photometer, with the ability to ob- 
serve at 75 /im, 110 /im, or 175 /im, will provide a measure 
of Hll region SEDs at wavelengths near their peak emissions. 
The SPIRE camera will observe at 250 /im, 350 /im, and 
500 /im (albeit with a lower resolution than PACS) and can 
provide a better constraint on the Rayleigh- Jeans side of the 
dust blackbody emission. NASA's 2.5 m SOPHIA telescope 
will also have a suite of high-resolution instruments capable 



of observing star forming regions throughout the IR, includ- 
ing near the expected peaks of the Hll region SEDs. 

This work is based on observations made with the Spitzer 
Space Telescope, which is operated by the Jet Propulsion Lab- 
oratory, California Institute of Technology, under NASA con- 
tract 1407. The SAGE-LMC project has been supported by 
NASA/Spitzer grant 1275598 and Meixner's efforts have had 
additional support from NASA NAG5- 12595. 
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